Novel Mechanism of Supersolid of Ultracold Polar Molecules in Optical Lattices 
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We study the checkerboard supersolid of the hard-core Bose-Hubbard model with the dipole- 
dipole interaction. This supersolid is different from all other supersolids found in lattice models 
in the sense that superflow paths through which interstitials or vacancies can hop freely are absent 
in the crystal. By focusing on repulsive interactions between interstitials, we reveal that the long- 
range tail of the dipole-dipole interaction have the role of increasing the energy cost of domain wall 
formations. This effect produces the supersolid by the second-order hopping process of defects. We 
also perform exact quantum Monte Carlo simulations and observe a novel double peak structure 
in the momentum distribution of bosons, which is a clear evidence for supersolid. This can be 
measured by the time-of-flight experiment in optical lattice systems. 

PACS numbers: 03.75.Hh, 05.30.Jp, 67.85.-d 



Recently, there are great experimental efforts to cre- 
ate ultracold dipolar molecules in optical lattices (Tj-Q- 
In these systems, owing to the dipole-dipole interaction, 
various quantum phases are predicted theoretically, such 
as supersolid, devils staircase, crystals with long period- 
icities and so on 4~8j]. Experimentally, the realization of 
supersolids has attracted great interest. In recent study 
by quantum Monte Carlo calculations Q, several charac- 
teristic supersolid states, such as checkerboard-type and 
star-type supersolids, were found in hardcore bosonic sys- 
tems with the dipole-dipole interaction on a square lat- 
tice. Since optical lattice is one of ideal tools that can 
be compared with quantum Monte Carlo results directly, 
this system is a promising candidate for the realization 
of supersolids. 

Although the detection of the checkerboard supersolid 
in hard-core bosonic systems is highly plausible in ultra- 
cold polar molecules, its mechanism is not clear. The 
possibility of supersolid state itself has been already dis- 
cussed by quantum Monte Carlo calculations for bosonic 
systems on several lattices such as triangular lattice [a- 
[ig |. face centered lattice [111], and square lattice p^4l4| . 
The mechanism for all these supersolids can be explained 
by the existence of superflow paths through which vacan- 
cies or interstitials can hop freely in crystals, thus Bose- 
condensing. In contrast to them, the checkerboard super- 
solid apparently has no such superflow paths in the hard- 
core boson case, because the particle hopping is restricted 
to the nearest-neighbor sites. In fact, it has been estab- 
lished that no supersolid phase appears when the sys- 
tem has only the nearest-neighbor and the next-nearest- 
neighbor repulsions on the square lattice Therefore, 
the long -range tail of the dipole-dipole interaction seems 
to play an important role to stabilize the supersolid state. 
However, its role has not been fully discussed yet. 

In this paper, we study the mechanism of the checker- 
board supersolid in dipolar hard-core bosons in 2D opti- 
cal lattices. By discussing the stability of the supersolid 



against the domain wall formation, we reveal that the 
long-range tail has the effect of increasing the energy cost 
of the domain wall formations and change the mechanism 
of the supersolid. We also perform the exact quantum 
Monte Carlo simulations to investigate the properties of 
this novel supersolid. As a result, we observe the two 
types of peaks in the momentum distribution of bosons, 
which indicate superfluidity and solidity respectively. 

We consider the hardcore Bose-Hubbard model with 
the dipole-dipole interaction on a square lattice. This 
system can be realized by bosonic polar molecules in a 
2D optical lattice. In this paper, electric dipole-dipole 
moments are assumed to be arranged perpendicular to 
the 2D plane by applying the electric field. Thus the 
dipole-dipole interaction is expressed as 1/r 3 here. The 
Hamiltonian is given by 



H = -tY J {b%+h.c.) + V 

(id) 
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where bi(b\) is the annihilation(creation) operator of 
hardcore bosons on the site i, rii = b\bi is a particle 
number operator, t, V, and /i indicate the hopping pa- 
rameter, the strength of the dipole-dipole interaction, 
and the chemical potential respectively. r,-j is the dis- 
tance between the site i and j. In what follows, we take 
the lattice spacing as the unit of distance. means 
nearest-neighbor pairs. The system size is defined by 
N = L x L and the periodic boundary condition is im- 
posed. In this model, the existence of the checkerboard 
supersolid phase is shown by the quantum Monte Carlo 
simulations Q. 

Before we discuss the mechanism of the checkerboard 
supersolid in the model (TTJ), we explain that this super- 
solid is different from previous ones which have been 
found in the other lattice models by exact quantum 
Monte Carlo simulations [8HT^| . All previous supersolids 
are related to the presence of superflow paths in crys- 
tals. Here, we define superflow paths as paths on which 
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FIG. 1: (Color online) Examples of superflow paths for inter- 
stitials in crystals. The red circles and the orange circles rep- 
resent the bosons forming the crystal and interstitial bosons 
respectively. The paths are denoted by thick lines. The ar- 
rows indicate the hopping process of interstitials. (a) The 
filling factor p = 1/3 crystal of hard-core bosons with the 
nearest-neighbor repulsion on a triangular lattice, (b) The 
checkerboard crystal of soft-core bosons with the on-site re- 
pulsion U and the nearest-neighbor repulsion V nn on a square 
lattice. For simplicity, the system is depicted along the one- 
dimensional direction. The wavy lines represent the repul- 
sive interactions between the interstitial and the checkerboard 
background. 



defects such as interstitials or holes in crystals can hop 
strictly or approximately by the first-order perturbation 
of t. Through these paths, defects in crystal can hop 
freely and Bose-condense. In Fig. [TJ we show the exam- 
ples of superflow paths. If we consider hard-core bosons 
on a triangular lattice with the nearest-neighbor repul- 
sion, the bosons form a superlattice like Fig. [T] (a) at the 
1/3-filling In this case, the superflow path is the 

honeycomb lattice formed by the vacant sites. As we can 
see in this example, sublattices have the separate roles of 
forming the crystal and superfluidity in hard-core bosonic 
systems. Another example is soft-core bosons on a square 
lattice with the nearest- neighbor repulsion V nn , where the 
checkerboard crystal appears like Fig. [T] (b) at the 1/2- 
filling [13]. In this case, the superflow path is formed all 
over the lattice for interstitials if the on-site interaction U 
is nearly equal to zV nn , where z is the coordination num- 
ber. In other words, the paths allow the interstitials to 
hop on the checkerboard backgroud. The superflow paths 
are blocked in the hard-core limit U — > oo. However, the 
checkerboard supersolid of hard-core bosons can be sta- 
bilized in the presence of the dipole-dipole interaction as 
mentioned above. Therefore, the present supersolid is 
qualitatively different from previous ones and expected 
to have some other mechanism. 

We now discuss the checkerboard supersolid in our 
model. In Fig. [3J we show the checkerboard crys- 
tal doped with L/2 bosons. As illustrated in Fig. 
[H (a), supersolid can be realized by delocalization of 
defects [Hf, [l6[. If the system has only the nearest- 
neighbor interaction, domain wall formations like Fig. [5] 
(b) are energetically more favorable than the supersolid 
for small hopping amplitudes, as discussed in Ref. [l3| . 
This is due to two important features in short-range inter- 
acting systems. One is that the energetic cost for doping 
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FIG. 2: (Color online) The checkerboard crystal doped with 
interstitials. The red circles, the orange circles, and the ar- 
rows have the same meanings in Fig. [1] The dashed line 
and the wavy line represents a domain wall and one of the 
repulsive interactions between interstitials respectively, (a) 
Supersolid where interstitials are delocalized. (b) Domain 
wall where interstitials are aligned in a particular direction. 
In long-range interacting systems, interstitials can gain the 
kinetic energy by the second-order hopping process in both 
cases. Furthermore, repulsive interactions between intersti- 
tials appear. 



a single boson is independent of positions of the empty 
sites. The other is that interstitials forming domain walls 
can gain the kinetic energy by the first-order hopping pro- 
cess, whereas bosons can delocalize by the second-order 
hopping process. By these features, with doping of the 
small particle density ~ l/£, the supersolid state(Fig. [5] 
(a)) and the domain wall state(Fig. H](b)) have the same 
energy in the zeroth-order in i, whereas the domain wall 
state has a lower energy in the first-order in t. Conse- 
quently, the supersolid becomes unstable to the domain 
wall formation. 

However, the situation should be changed when the 
system has the long-range interactions, because the re- 
pulsive interaction does work between interstitials. As 
we can see in Fig. [5] (b) , interstitials forming the domain 
wall feel repulsive energies with each other in the pres- 
ence of the third or higher-nearest neighbor repulsion. By 
these repulsions, it costs higher energy than the config- 
uration where interstitials randomly occupy empty sites 
like Fig. [5] (a) even in the zeroth-order of t. In the do- 
main wall case, the kinetic energy gain derived from the 
second-order hopping process apparently exits. Since the 
gain is the same order in t as that in the case shown in 
Fig. H](a), it is naively expected that configuration en- 
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ergy of interstitials decides everything; the system prefers 
the supersolid state energetically. 

In order to confirm quantitatively that the long-range 
interactions are relevant and change the stability of the 
supersolid against the domain wall, we compare the en- 
ergies of Fig. [5] (a) and (b) by the perturbative calcu- 
lations. For simplicity, we consider the dipole-dipole in- 
teraction up to the fourth-nearest neighbor interactions. 
Adding a single boson to the checkerboard crystal costs 
the energy E n = zV + 2zV/(^/E) 3 — fi. For the domain 
wall configuration (b), by the repulsive interactions be- 
tween interstitials and the shift of the lower half configu- 
ration, the energy is raised by 2V/(\/5) 3 — V/2 3 ~ 0.054V 
per interstitial. When we introduce the kinetic terms, 
we obtain the kinetic energy as ~ — 14i 2 /V by the the 
second-order perturbative calculation. By these terms, 
the energy cost for forming the domain wall state is es- 
timated as E D w ~ E + 0.0541^ - Ut 2 /V per intersti- 
tial. As shown in Ref. [7|, the checkerboard supersolid 
phase is observed around the checkerboard Mott lobes for 
t/V > 0.05. When we use the small hopping parameter 
t/V = 0.05 for estimation of the energy above, the zeroth 
energy 0.0541/ is still larger than the kinetic energy gain 
Ut 2 /V ~ 0.035V. Thus, we obtain that E DW is larger 
than Eq; Edw > Eq. This is a result of higher configura- 
tion energy and lower kinetic energy gain than the case 
with only the nearest-neighbor repulsion. On the other 
hand, for the supersolid state (a), the repulsive energies 
between the interstitials are negligible, because the inter- 
stitials are very far away from each other in the thermo- 
dynamic limit. Hence, the zeroth term is only Eq. In the 
presence of the kinetic term, the energy is lowered by the 
second or higher-order hopping process. Thus, the energy 
cost of a single interstitial is Ess = Eq — 0(t 2 /V), which 
is smaller than Eq; Ess < Eq. From these estimations, 
we conclude that the supersolid becomes stable against 
the domain wall formation owing to the effect of the re- 
pulsive interactions between interstitials. The above esti- 
mation does not change up to the presence of the eighth- 
nearest-neighbor repulsion, because the second-nearest 
interstitials interact by the ninth-nearest-neighbor repul- 
sion which is negligibly small in the estimation above. 
More specifically, the third and fourth-nearest-neighbor 
interactions play the dominant role of increasing the en- 
ergy cost of the domain wall formation. Because of the 
absence of these interactions, the checkerboard supersolid 
have not been found in systems with only the nearest 
neighbor repulsions and the next-nearest neighbor repul- 
sions. Although we have discussed the interstitial-based 
supersolid above, the vacancy-based supersolid is also 
possible by the same mechanism because of the particle- 
hole symmetry in hard-core bosonic systems. 

We next perform exact quantum Monte Carlo simu- 
lations to investigate the nature of the supersolid phase 
described above. The notorious difficulty in simulations 
of long-range interacting systems is the computational 
cost 0(N 2 ), which is O(N) if the system has only short- 
range interactions. In order to overcome this difficulty, 
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FIG. 3: (Color online) Off-diagonal correlation function 
Csf(t") and diagonal correlation function Ccdw(?") as a func- 
tion of the distance along the a;-axis. CsF(r) and CcDw(f ) are 
shown as logarithmic plots and semi-logarithmic plots respec- 
tively. Error bars are drawn, but most of them are smaller 
than the symbol size. 
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FIG. 4: (Color online) Momentum distributions of bosons 
nsF(fc x ,fcy) at L = 120. The other parameters are the same 
as in Fig. [3] The central sharp peak is located at fc = 0. The 
four smaller peaks are at k — (ji, tv) and its equivalent points 
in the momentum space. 

we use a new hybrid algorithmf!?) of modified directed 
loop algorithm fisj and 0(N) Monte Carlo method (l9j. 
In our simulations, we applied the Ewald summation 
method[20| to avoid any cutoff dependence of the long- 
range interactions. 

In Fig. [3j we show the off-diagonal correlation func- 
tion CsF{fij) = (bibj) and the diagonal correlation func- 
tion CcDw( r ij) = ( n i n j) — ( n i) 2 as a function of the 
distance along the x-axis. The parameters are chosen 
as (t/V,fj,/V,T/t) = (0.15,3.3,0.1) where the vacancy- 
based checkerboard supersolid phase appears. To ob- 
serve the long-distance behavior of the correlation func- 
tions, we perform simulations up to the system size 
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120 x 120, which needs high computational cost with- 
out the 0(N) method. From this result, we can find 
that Csf (r) shows a power-law decay, implying the pres- 
ence of an off-diagonal quasi-long-range order. Simulta- 
neously, Ccdw(^) shows the presence of a diagonal long- 
range order. Thus, this phase is the supersolid phase 
in a broad sense. As the most characteristic feature in 
this supersolid, we have observed a clear zig-zag pattern 
of the power-law decay in Csf(t)- Obviously, this pat- 
tern is a sign that the superfluid component is affected 
by the checkerboard-type structure. In Fig. |H we show 
the Fourier transformation of the off-diagonal correlation 
function n SF (k) = 1/YV^- C SF {r ij )e lk - r ^ . As in the 
conventional superfluid, we can observe the sharp peak 
at k = 0, which means superfluidity. In addition to this, 
as a result of the zig-zag pattern, we can also observe 
the other peaks at k = (tt, tt) and its equivalent points, 
although it is much weaker than the former one. The cor- 
responding peaks can be measured by the time-of-flight 
experiment and the observation of the two types of peaks 



will be a clear evidence for supersolidity. 

In summary, we have revealed the novel mechanism of 
the supersolid of hardcore bosons with the dipole-dipole 
interaction. With repulsive long-range interactions, the 
superflow paths are not required any more for realizing 
supersolids, and the mechanism by the second-order hop- 
ping process of interstitials becomes relevant. Since the 
checkerboard supersolid is due to the second-order hop- 
ping process, the transition temperature T c should be 
scaled as T c oc i 2 /V , and it can be very low in real ex- 
perimental settings. 
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